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Quantum and Wave Dynamical Chaos in Superconducting Microwave Billiards 
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Experiments with superconducting microwave cavities have been performed in our laboratory for 
more than two decades. The purpose of the present article is to recapitulate some of the highlights 
achieved. We briefly review (i) results obtained with flat, cylindrical microwave resonators, so-called 
microwave billiards, concerning the universal fluctuation properties of the eigenvalues of classically 
chaotic systems with no, a threefold and a broken symmetry; (ii) summarize our findings concerning 
the wave-dynamical chaos in three-dimensional microwave cavities; (iii) present a new approach for 
the understanding of the phenomenon of dynamical tunneling which was developed on the basis of 
experiments that were performed recently with unprecedented precision, and finally, (iv) give an 
insight into an ongoing project, where we investigate universal properties of (artificial) graphene 
with superconducting microwave photonic crystals that are enclosed in a microwave resonator, i.e., 
so-called Dirac billiards. 

PACS numbers: 03.65.Sq, 05.45.Mt, 11.30.Er 41.20.Jb, 73.40.Gk, 74.50.tr 


Quantum chaos, a field that focuses on the 
quantum manifestations of classical chaos, at¬ 
tracted a lot of attention during the last 35 
years. The signatures of classical chaos in uni¬ 
versal properties of the corresponding quantum 
system are indeed of relevance in nuclear physics, 
atomic physics, condensed matter physics, op¬ 
tics and in acoustics. They were observed in the 
statistical properties of the eigenvalues and the 
wave functions, in the fluctuation properties of 
the scattering matrix elements of chaotic scat¬ 
tering processes, in the transport properties of 
quantum dots and also in systems, where time- 
reversal invariance is broken, e.g., by a magnetic 
field. Billiards provide a most appropriate sys¬ 
tem for the study of features from the field of 
quantum chaos, because the degree of chaoticity 
of their classical dynamics only depends on their 
shape. The eigenvalues and wave functions of 
two-dimensional quantum billiards are measured 
since two decades in our laboratory using flat 
cylindrical, superconducting and normal conduct¬ 
ing microwave resonators, respectively. Thereby 
large sequences of approximately 1000 eigenval¬ 
ues were determined in the first experiments with 
a hitherto never achieved high precision yield¬ 
ing new insight into quantum chaos phenomena, 
because at that time their computation was ei¬ 
ther not possible or only feasible with a large 
numerical effort. We will present the results of 
these and subsequent experiments and also con¬ 
cerning our latest studies of dynamical tunnel¬ 
ing. The latter demonstrate that even nowa¬ 
days the measurements with superconducting mi¬ 
crowave resonators may provide valuable new in- 
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formation that is indispensable for the develop¬ 
ment of a theoretical description of quantum phe¬ 
nomena like dynamical tunneling, even though 
eigenvalues may now be computed with high ac¬ 
curacy. Most of the works on quantum chaos con¬ 
centrate on quantum systems described by the 
non-relativistic Schrodinger equation. Only a few 
studies focus on the properties of the eigenval¬ 
ues and the wave functions of relativistic billiards 
that are governed by the Dirac equation of a 
massless spin-1/2 particle. We started lately with 
their experimental investigation and will present 
the system and recent results on the statistical 
properties of the associated eigenvalues in this ar¬ 
ticle. 


I. INTRODUCTION 

This article provides a brief survey on groundbreak¬ 
ing experiments that were performed in the last two 
decades in our group with superconducting microwave 
billiards. They focused on the understanding of the 
generic properties of quantum systems like nuclei and of 
condensed matter structures like graphene or fullerene. 
In the first case, we exploited the equivalence of the 
scalar Helmholtz equation governing flat, metallic mi¬ 
crowave resonators, so-called microwave billiards, and the 
non-relativistic Schrodinger equation of the correspond¬ 
ing quantum billiard, to address problems from the fields 
of quantum chaos and compound nucleus reaction the¬ 
ory m- In the second case, the analogy between the 
scalar Helmholtz equation describing microwave photonic 
crystals squeezed between two metal plates or enclosed 
in a microwave billiard, so-called Dirac billiards, and the 
Dirac equation was used to experimentally investigate 
relativistic phenomena occurring in graphene. Some de¬ 
tails on the experimental setups for the measurements of 
resonance spectra at superconducting conditions will be 
provided in Sec. HIl 
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FIG. 1: Part of the transmission spectrum, i.e., the outcou- 
pled over incoupled power Pout/Pin versus the excitation fre¬ 
quency, of a microwave billiard with the shape of a quarter 
stadium measured at room temperature (upper panel) and 
at 2 K (lower panel), where it is normal conducting and 
superconducting, respectively. The shape of the resonator 
and the positions of the antennas are illustrated in the inset. 
Reprinted from Phys. Rev. Lett. 69, 1296 (1992). 


A classical billiard [3] is a bounded domain, in which a 
pointlike particle moves freely and is reflected elastically 
on impact with the boundary. Since its dynamics only 
depends on its shape, such systems are widely used to 
investigate the effects of the classical dynamics on prop¬ 
erties of the associated quantum system me]- The cor¬ 
respondence between the latter, i.e., a quantum billiard 
and a microwave billiard of the same shape, holds be¬ 
low a certain excitation frequency / of the microwaves 
that are sent into the resonator, / < fmax = c/2d^ 
where the electric field vector is parallel to the cylinder 
axis. Here, c is the velocity of light and d is the height 
of the resonator. Originally, this analogy was used to 
study universal fluctuation properties of the eigenvalues 
of a quantum billiard with microwave billiards at room 
temperature IHHII]. The resolution of the resonances in 
the spectra, however, was very limited due to the rela¬ 
tively low quality factor Q 10^ for normal conduct¬ 
ing microwave billiards. We obtained the eigenvalues of 
(non) relativistic quantum billiards by measuring reso¬ 
nance spectra to a very high accuracy with supercon¬ 
ducting microwave billiards with quality factors of up to 
10^. This is an indispensable requirement for the deter¬ 
mination of complete sequences of several hundreds and 
even thousands of eigenvalues, as illustrated in Fig. 
Shown are transmission spectra of a microwave billiard 
made out of niobium metal with the shape of a quar¬ 
ter stadium (see Sec. Im]) measured at room temperature 
(upper panel), where the resonator is normal conducting 
and at 2 K, where it is superconducting. Its eigenfre- 
quencies, and thus the eigenvalues of the corresponding 
quantum billiard, are obtained from the positions of the 


resonances. In the upper spectrum in Fig. [^many reso¬ 
nances overlap, whereas all resonances are well resolved 
in the lower one, so that the determination of the com¬ 
plete sequence of eigenvalues is achievable only in the su¬ 
perconducting case. Because of this capability we were 
able to tackle a wide range of problems that are far from 
feasible with normal conducting microwave billiards. 

Since the resonance spectra are measured by emitting 
microwave power into the resonator via one antenna and 
receiving it at the same or another one, microwave bil¬ 
liards can also be viewed as scattering systems. We uti¬ 
lized this capability and studied the universal properties 
of quantum scattering systems with an unprecedented 
accuracy. Due to limitation in space we focus here only 
on a few of our experiments. In Secs. [ITT] and [V^ we 
present results on the spectral properties of classically 
chaotic systems. In Sec. IIVI we review our studies con¬ 
cerning the effect of symmetry breaking on the proper¬ 
ties of the eigenvalues and the wave functions of chaotic 
systems. We also investigated the spectral properties of 
three-dimensional resonators that do not have a quantum 
analogue. These experiments provided additional impor¬ 
tant information to previously performed ones [IMS] 
of, due to the superconductivity, hitherto never achieved 
precision. One example will be discussed in Sec. [Vl 

The analysis of the spectral properties addresses the 
manifestation of the gross structure of the Poincare sur¬ 
face of section (PSOS), namely the relative sizes of the 
regions of regular and chaotic dynamics, in the associated 
quantum system. The phenomenon of quantum dynami¬ 
cal tunneling [T6|, that is, the classically forbidden tran¬ 
sition between states localized in separated parts of the 
classical phase space has been studied theoretically inten¬ 
sively, but only in a few experiments [MU], because the 
measurements require an exceptionally high precision. In 
Sec. Eni we will review the main results of our investi¬ 
gations of chaos-assisted tunneling m and of dynamical 
tunneling m- 

Only lately we started experimental investigations of 
relativistic phenomena associated with graphene and will 
present some results of our most recent experiments 

in Sec. [ylTIl 


II. SOME DETAILS OF THE EXPERIMENTAL 
SETUPS 

As emphasized in Sec. [Tjand illustrated in Fig.[^for the 
determination of complete sequences of eigenvalues of a 
quantum billiard from the resonance spectra of the cor¬ 
responding microwave billiard, high-precision measure¬ 
ments are indispensable. Therefore, we used microwave 
billiards that were superconducting at the temperature 
of liquid helium T = 4.2 K. The first billiards were con¬ 
structed in the CERN workshop in Geneva. They were 
made from niobium and quality factors Q ~ 10^ — 10^ 
were achieved. In the very first experiments a microwave 
billiard with the shape of a quarter stadium, reported on 
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in Sec. Eni was inserted into one of the cryostats of the 
superconducting Darmstadt electron linear accelerator S- 
D ALIN AC [20] and the measurements were performed at 
2 K. A part of a spectrum is shown in the lower panel 
of Fig. 12 In the subsequent experiments the microwave 
billiards were cooled down to 4.2 K with liquid helium 
in a separate LHe bath cryostat. The niobium cavities 
were made from one piece, so it was not possible, e.g., to 
vary their shapes or to introduce a scatterer into them. 
Therefore, since about the year 2000 the resonators were 
constructed from three copper parts, a bottom plate, a 
middle one with a hole of the shape of the billiard and a 
lid. Recently, we have constructed them from two brass 
plates, where a bassin in the shape of the billiard was 
milled out of the bottom one. To attain superconductiv¬ 
ity all parts of the resonator were lead plated and then 
squeezed on top of each other tigthly with screws along 
the boundary of the billiard. 

For the measurement of the resonance spectra mi¬ 
crowave power was emitted into the resonator via one 
wire antenna and coupled out either via the same an¬ 
tenna or via a second one, thus yielding a reflection and a 
transmission spectrum, respectively. A vectorial network 
analyzer determined the relative phase and amplitude of 
the output and input signal. This provided the scatter¬ 
ing matrix elements S^a for fhe scattering from antenna 
a to antena 6, with a,b = 1, 2, • • • depending on the num¬ 
ber of antennas introduced into the resonator via holes in 
the top plates. As already mentioned above, the height 
d of the resonator determines the maximum frequency 
fmax = c/2d below which the analogy to the quantum 
billiard of corresponding shape holds. The billiard’s size, 
and thus the number of attainable eigenmodes [21], is 
bounded by that of the cryostat. Accordingly, in the 
first measurements with microwave billiards that had a 
height d = 7 mm approximately 1000 eigenvalues could 
be determined below fmax — 20 GHz, the most recent 
measurements yielded several thousand eigenvalues using 
resonators with a height of d = 3 mm which corresponds 
to a maximum frequency fmax — 50 GHz. In sections HI, 
VI, Vllamd VHI we will present results from all types of 
resonators. 


III. QUANTUM CHAOS IN CLOSED AND 
OPEN SYSTEMS: THE BUNIMOVICH STADIUM 
BILLIARD 


A multitude of the experiments addressed the statisti¬ 
cal properties of the eigenvalues and the wave functions 
of billiards of different shapes. They focused on problems 
related to the field of quantum chaos, as, e.g. the validity 
of the Bohigas-Giannoni-Schmit conjecture [llllEllISS] 
for chaotic systems. It states that the spectral fluctua¬ 
tion properties of a chaotic system coincide with those of 
random matrices from the Gaussian orthogonal ensemble 
(GOE), if time-reversal (T) is conserved, and with those 
of random matrices from the Gaussian unitary ensem¬ 


ble (GUE), if it is violated. If the classical dynamics is 
regular they follow the Poisson statistics. There also ex¬ 
ist random-matrix theory (RMT) predictions for systems 
with mixed regular-chaotic dynamics, and for chaotic 
systems with a partially broken symmetry, as, e.g., the 
isospin mixing observed in nuclear physics (see Sec. Ej), 
and with partially violated Tinvariance. 

The first experiments were performed with a niobium 
microwave billiard with the shape of a quarter Buni¬ 
movich stadium [T] shown in Eig. The radius of the 
quarter circle equaled r = 200 mm and the length of the 
rectangular part a = 560 mm. 



FIG. 2: (Color online) Photograph of the two-dimensional 
microwave billiard with the shape of a quarter stadium. The 
resonator was constructed from niobium which is supercon¬ 
ducting below Tc = 9.2 K. Its height is 8 mm so the system 
is governed by the scalar Helmholtz equation below the fre¬ 
quency fmax ~ 18.9 GHz. The microwave billiard was cooled 
down to 2 K in the superconducting Darmstadt electron linear 
accelerator S-DALINAC and quality factors of Q ~ 10^ — 10^ 
were achieved. 


A complete sequence of 1060 resonance frequencies was 
identified. For the study of the spectral properties the 
spectra were unfolded to mean resonance spacing unity 
with Weyl’s formula [21] which provides an analytical ex¬ 
pression for the smooth part of the integrated resonance 
density N{f) in terms of a polynomial of second order, 
Nsmoothif) = - Bl{4c)f + C, with A the 

area and B the perimeter of the billiard. Figure]^ shows 
the fluctuating part of the integrated resonance density, 
Nfiucif) = N{f)-Nsmooth{f) which exhibits fast fluctu¬ 
ations and smooth oscillations. The latter are due to non¬ 
generic periodic orbits (POs), so-called bouncing-ball or¬ 
bits (BBOs), that bounce back and forth between the two 
parallel side walls of the rectangular part of the stadium 
billiard. The red (gray) full line shows the analytical 
result for their contribution. 


NM = 





( 1 ) 

It was derived on the basis of a semiclassical method [24] . 
While the bouncing-ball orbits are of measure zero in the 
classical dynamics they have a strong effect on the spec¬ 
tral propertes of the quantum billiard, as demonstrated 
in Fig. Shown is the Dyson-Mehta As{L) statis¬ 
tics m, which gives the least-square deviation of the 
integrated density of the unfolded resonance frequencies 
from the best straight-line fit in the interval of length L. 
Since the classical dynamics of a Bunimovich stadium is 
completely chaotic [26] it is expected to coincide with the 
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RMT results for the GOE. The circles and the diamonds 
show the experimental results obtained when unfolding 
the resonance frequencies / by replacing them as usually 
by / = Nsmoothif), and by / = Nsmoothif) + Nbb{f) 
taking into account the effect of the BBOs, respectively. 
It is clearly visible that the agreement between the ex¬ 
perimental As statistics and the GOE is very good only 
in the latter case. 



FIG. 3: (Color online) Fluctuating part of the integrated 
density of the resonance frequencies obtained for the mi¬ 
crowave billiard in Fig. (black line) and the contributions 
of the bouncing ball orbits (red line) deduced from Eq. 0. 
Reprinted from Phys. Rev. Lett. 69, 1296 (1992). 



FIG. 4: The As statistics of the resonance frequencies of the 
microwave billiard in Fig. [^before (diamonds) and after (cir¬ 
cles) extraction of the contribution from the bouncing-ball 
orbits compared to the RMT result for the GOE (full line). 
Reprinted from Phys. Rev. Lett. 69, 1296 (1992). 

As was noted above, for the measurement of the reso¬ 
nance spectra microwave power is emitted into the res¬ 
onator via one antenna thereby exciting an electric field 
mode in its interior and an output signal is received at 
the same or another antenna. Hence, the resonator is an 
open scattering system with the antennas acting as single 
scattering channels. The associated S matrix comprises 
the scattering between the input and the output anten¬ 
nas and dissipation into the walls, which is accounted for 
by introducing a large number of weakly coupled ficti¬ 
tious scattering channels. The ^'-matrix formalism for 


compound nucleus reactions m is equivalent to that 
for microwave resonators [28]. We used this analogy to 
investigate universal properties of the associated S ma¬ 
trix and more generally for quantum scattering processes 
with intrinsic chaotic dynamics. In regions of isolated 
resonances, as is the case in the measurements with the 
superconducting quarter-stadium billiard, the ^'-matrix 
element for the scattering from antenna a to antenna b 
is given as 

Sbaif)=Sba-ij^^^^^ ( 2 ) 

in the vicinity of the resonance ji. Here, V and F^^ are 
the partial widths related to the emitting and the receiv¬ 
ing antennas a and 6, respectively, and F^ is the total 
width of the resonance, F^ = Yfc=i It in¬ 

cludes the partial widths of the 3 antennas attached to 
the microwave billiard (see Fig.|^ and the dissipation in 
the walls. Note, that only for the spectra that were mea¬ 
sured at 2 K, i.e., for those of the quarter-stadium billiard 
all resonance parameters could be determined from a fit 
of the Lorentzian in Eq. to the resonances [29]. The 
reason is, that the contributions T^^^iss from the losses 
due to the dissipation in the walls were negligible as com¬ 
pared to those from the partial widths, whereas this is 
no longer the case in the measurements at liquid Helium 
temperatures of 4.2 K, see below. Eigurej^ shows the to¬ 
tal widths F^ and the partial widths F ^2 associated with 
antenna 2 versus the resonance frequencies /^. The data 
fluctuate strongly around a slow secular variation which 
is well approximated by a polynomial of 5th order and 
was removed for the analysis of their statistical properties 
by scaling the data with the latter. 

The distribution of the partial widths is expected to 
coincide with that of a xt distribution with the number 
of degrees of freedom = 1, which is a Porter-Thomas 
distribution. The sum of 2 and 3 partial widths should 
coincide with that of 2 and 3 degrees of freedom, respec¬ 
tively, if the partial widths indeed are uncorrelated. As 
demonstrated in Fig. the agreement between the ex¬ 
perimental and the expected results is very good [61130]. 

Furthermore, the ^'-matrix autocorrelation function 
was computed using the experimental resonance spectra 
and the RMT predictions for those of isolated resonan¬ 
ces [31] were confirmed. The Fourier tranforms of the 
^'-matrix autocorrelation function revealed an algebraic 
decay. This actually is expected for isolated and weakly 
overlapping resonances [30l|3l]. 

We would like to close this section with a few remarks 
concerning our experimental tests of analytical results 
within the theory of chaotic scattering which has been 
largely developed in the framework of the statistical the¬ 
ory of compound-nucleus reactions [2]. Predictions for 
the fluctuation properties of the associated S matrix 
were obtained based on RMT and on the supersymmetry 
method [32ll35] . To be more precise, analytic expressions 
for ^'-matrix correlation functions and the distribution of 
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a microwave resonator with the shape of a chaotic tilted 
stadium billiard allowing tests also with unprecedented 
accuracy in the regions of isolated and weakly overlap¬ 
ping resonances [34l [36H39] . Experiments were further¬ 
more performed with microwave billiards where a partial 
Tviolation was induced by a magnetized ferrite [40]. For 
the analysis of the latter we extended predictions derived 
for T-invariant systems by Verbaarschot, Weidenmiiller 
and Zirnbauer [32] for the ^'-matrix autocorrelation func¬ 
tions to systems with partial Tviolation. Due to the large 
sets of scattering data measured as function of the exci¬ 
tation frequency we could also study higher order corre¬ 
lation functions, as e.g. the cross-section autocorrelation 
function to test RMT predictions on cross-section fluctu¬ 
ations [371439] . 


FIG. 5: The total widths (upper panel) and the partial widths 
(lower panel) of the resonances measured with the super¬ 
conducting microwave billiard shown in Fig. [^ versus their 
frequencies. The full lines are polynomial fits to the data, 
that describe the secular dependence of the widths on the 
frequency. Reprinted from Phys. Rev. Lett. 74, 62 (1995). 



FIG. 6: Distributions of the partial widths F^a associated 
with antennas a = 1, 2, 3 normalized by their respective mean. 
Those of F^i (upper panel) and of the sum of F^i + F ^2 
(middle panel) and of all (lower panel) partial widths (his¬ 
tograms) of the resonances measured with the superconduc¬ 
ting microwave billiard in Fig. [^ respectively, are displayed. 
The full curves show the xl distributions of z/ = 1, 2 and 3 de¬ 
grees of freedom from top to bottom. Adopted from Ref. [3. 
With kind permission from Springer Science and Business Me¬ 
dia. 

the S'-matrix elements have been derived which are valid 
from the regime of isolated resonances to that of strongly 
overlapping ones. For their experimental verification, 
scattering data where taken at room temperature using 


IV. STUDY OF SYMMETRY BREAKING 
WITH COUPLED QUARTER-STADIUM 
BILLIARDS 

Conservation laws, that is, the invariance of physical 
quantities with respect to certain transformations, play 
an important role in nearly all fields of physics. In 1918, 
Emmy Noether formulated a theorem [41], which states 
that any differentiable symmetry of the action of a phys¬ 
ical system has a corresponding conservation law. In 
a quantum system each conserved quantity corresponds 
to a quantum number. We investigated the symmetry 
breaking in a system that is characterized by two quan¬ 
tum numbers, like, e.g., that of the isospin symmetry in 
^^Al or which is broken due to the Coulomb inter¬ 
action, see [42]444] and references cited therein. In ^^Al 
a complete sequence of 107 states of positive parity from 
the ground state up to about 8 MeV excitation energy 
was identified. Among these 75 had isospin quantum 
number T = 0 and 32 T = 1, so one has the situation of 
two symmetry classes with each characterized by a GOE, 
since the corresponding classical dynamics is chaotic. 
The upper panel of Fig. [^ shows their nearest-neighbor 
spacing distribution (histogram). It should coincide with 
that of two uncoupled GOEs, plotted as dashed line in 
Fig-IZl if the isospin were a good quantum number. If, on 
the other hand, the symmetry were completely broken, 
that is, if there were no isospin symmetry, the distribu¬ 
tion should coincide with that of one GOE, shown as 
dotted line. The fact, that the nearest-neighbor spac¬ 
ing distribution coincides with neither of the two clearly 
indicates that the isospin is partially broken in ^^Al, or 
equivalently, that the two sets of states with isospin quan¬ 
tum numbers T = 0 and T = 1 are mixed. This mixing 
of the states of two chaotic systems is modelled by an 
ensemble of random matrices of the form 



which is a special case of the Rosenzweig-Porter 
model [45]. Here, the first part preserves the symme- 
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FIG. 7: Nearest-neighbor spacing distribution of the energies 
of the first 107 excited states with positive parity of A1 (up¬ 
per panel) and of the resonance frequencies of the uncoupled 
stadium billiards (lower panel). The dashed lines correpond 
to the curve for two uncoupled GOEs, the dotted ones to that 
of one GOE. 

try, where, in the case of ^^Al, Hi and H 2 are GOE ma¬ 
trices and describe the spectral properties of the states 
with T = 0 and T = 1, respectively. The symmetry 
breaking effect of the Coulomb interaction, that is, the 
mixing of the isospin states is accounted for by the sec¬ 
ond term, with D denoting the mean level spacing. For a 
non-vanishing symmetry-breaking parameter A it yields 
a coupling of Hi and H 2 . 

The symmetry breaking actually could be tested ex¬ 
perimentally by mapping this RMT model onto a sys¬ 
tem consisting of two coupled Bunimovich stadium bil¬ 
liards [46] shown schematically in Fig. A smaller mi- 


7i=l 72=1-0 



400mm 560mm 


FIG. 8: Schematic view of the two coupled Bunimovich sta¬ 
dium billiards. The locations of the antennas are marked by 
circles. Reprinted from Phys. Rev. Lett. 81, 4847 (1998). 

crowave stadium billiard with the ratio 71 = 1.0 between 
the length of the rectangle and the radius of the quarter 
circle was coupled to the one presented in the previous 
section III and shown in Fig.j^which had a ratio 72 = 1.8. 
Thus, the number of resonances identified up to 17.5 GHz 
in the former, Ni = 608, is smaller than that in the lat¬ 
ter, N 2 = 883. The dynamics in the two uncoupled cav¬ 


ities is simulated in the RMT model [Eq. by GOE 
matrices Hi and H 2 with the dimensions equal to Ni 
and 5 respectively. For the coupling, the microwave 
billiards were put on top of each other and a niobium 
pin, 2 mm in diameter, penetrated through holes in their 
walls into both of them. Its size was controlled by the 
associated penetration length. This led to the formation 
of a TEM mode between the walls and the pin. The mix¬ 
ing of the Ni and N 2 modes of both microwave billiards 
via the TEM mode was modelled by an Ni dimensional 
vector V and an N 2 dimensional vector w with Gaussian- 
distributed random numbers as entries, respect iv^, thus 
yielding a dyadic structure for V = vw^ in Eq. The 
distribution of its matrix elements is given by a modified 
Kq Bessel function. For a more detailed description see 
Refs. j46l[ 47 ] . 

Altogether 6 transmission spectra were measured with 
three antennas for 5 different couplings. The effect of 
the coupling on the spectra is illustrated in Fig. The 
change of the coupling led to a shift of resonances, whose 
associated electric field strength distribution, i.e., wave 
function, was non-vanishing at the position of the pin. 
We investigated the influence of the symmetry break¬ 
ing on the eigenvalues of a chaotic system in terms of 
the spectral fluctuation properties of the resonance fre¬ 
quencies, and that on the wave functions in terms of the 
resonance strengths r^^F^^ in Eq. (I^. The latter pro¬ 
vide information on the electric field strengths, i.e., the 
wave function components at the positions of the anten¬ 
nas. The lower panel of Fig.j^shows the nearest-neighbor 
spacing distribution of the unfolded resonance frequen¬ 
cies for the case of no coupling. It is well described by 
the RMT model [Eq. with the coupling parameter 
set to zero, A = 0 (dashed line). 



10.40 10.41 10.42 10.43 10.44 10.45 

f (GHz) 


EIG. 9: Part of the transmission spectra measured with no 
coupling, with weak coupling and with the strongest achieved 
coupling from top to bottom. The arrows mark a few reso¬ 
nances that are shifted due to the coupling. Reprinted from 
Phys. Rev. Lett. 81, 4847 (1998). 
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In order to determine the symmetry-breaking parame¬ 
ter for the different experimental realizations of the cou¬ 
pling strength we, firstly, fitted as outlined in Ref. [46] the 
analytical result [43l |48] for the number variance 
of the eigenvalues of the random matrices given in Eq. 
to that of the unfolded resonance frequencies. In Fig. 
the experimental data points with error bars and the an¬ 
alytical results (full lines) are shown for three different 
couplings together with those for the limiting cases of no 
symmetry, i.e., for 1 GOE (dashed line), and of symmetry 
conservation, i.e., 2 GOEs (dotted line). 




FIG. 11: Experimental strength distributions obtained from 
the resonances of the coupled stadium billiards for different 
couplings (histograms). The distributions of the analytical 
model derived from Eq. § including the threshold of detec¬ 
tion, best fitting the experimental results are shown as dashed 
lines. The symmetry breaking parameters deduced from the 
fits equal A = 0.003, 0.09, 0.14 and 0.20 from the upper left 
to the lower right panel. 

ing, r-^/D 0.25, i.e., approximately every fourth reso¬ 

nance is influenced by the coupling. 


V. WAVE-DYNAMICAL CHAOS IN 
THREE-DIMENSIONAL RESONATORS 


FIG. 10: The statistics of the measured resonance frequen¬ 
cies of the coupled stadium billiards for no coupling, a weak 
coupling and the strongest achieved coupling from top to bot¬ 
tom (filled circles). The dashed and the dotted lines show the 
RMT results for two uncoupled GOEs and one GOE, respec¬ 
tively. The full line results from a fit of the analytical result 
for the statistics of the model [Eq. §] to the experimental 
one. This yielded A < 0.03, A = 0.13 and A = 0.20 from top to 
bottom. Reprinted from Phys. Rev. Lett. 81, 4847 (1998). 


Secondly, in Ref. m, we compared the distributions 
of the resonance strengths obtained from the 

measured resonances with a qualitative model developed 
based on the eigenvectors of the random matrices in 
Eq. Figure 11 shows the experimental (histograms) 
and analytical (dashed lines) results for one antenna 
combination and for four different couplings. Here, we 
took into account the experimental threshold of detec¬ 
tion for narrow resonances and resonances with small 
strengths [49] . The symmetry-breaking parameters A ob¬ 
tained from both procedures agree within the experimen¬ 
tal errors. The symmetry breaking actually turned out 
to be relatively weak even in the case of the strongest 
coupling, but nevertheless to be of the same size as the 
isospin symmetry breaking in ^^Al, which yielded for the 
spreading width F'*^ with respect to the mean level spac- 


Another important issue besides the spectral fluctu¬ 
ation properties of a quantum system within the field 
of quantum chaos concerns the semiclassical approach. 
In this context, one of the most important achieve¬ 
ments was the periodic orbit theory (POT) pioneered by 
Gutzwiller m, which paved the way for an understand¬ 
ing of the impact of the classical dynamics on the prop¬ 
erties of the corresponding quantum system in terms of 
purely classical quantities. In contrast, RMT addresses 
exclusively system independent, i.e., universal proper¬ 
ties. Yet, only recently, the Bohigas-Giannoni-Schmit 
conjecture [5] was proven on the basis of Gutzwiller’s 
trace formula, which provides a semiclassical approxima¬ 
tion for the fluctuating part of the resonance density of 
a quantum system in terms of a sum over the POs of the 
corresponding classical system m- We have performed 
experimental tests of Gutzwiller’s trace formula and its 
analogue for systems of chaotic, mixed or regular dy¬ 
namics with metallic and dielectric resonators [IllSlHSS]. 
We also tested a Gutzwiller-like trace formula derived 
by Balian and Duplantier for three-dimensional electro¬ 
magnetic resonators [56ll58] . Here the main difficulty to 
overcome was the design of a microwave resonator that 
is completely chaotic and permits only a few non-generic 
modes. They were even present in our latest experiments 





















with the resonator shown in Fig. which had the shape 
of a three-dimensional stadium billiard. Note, that there 



FIG. 12: (Color online) Photograph of the resonator with the 
shape of a three-dimensional desymmetrized stadium billiard. 
It was constructed from niobium. 

is no analogy between three-dimensional microwave res¬ 
onators and quantum billiards. Still, the former are most 
suitable for the study of wave-dynamical phenomena in 
chaotic systems. 

The three-dimensional desymmetrized stadium bil¬ 
liard consisted of two quarter cylinders with radii ri = 
200.0 mm and r 2 = r\j\/2 mm that were rotated by 
90° with respect to each other so that the height of each 
cylinder coincided with the radius of the other one. The 
resonator was constructed from niobium and the mea¬ 
surements were performed at liquid helium temperature, 
yielding exceptionally high quality factors of up to 10^. 
Consequently, a complete sequence of 18764 resonance 
frequencies could be identified for excitation frequencies 
of up to 20 GHz. Parts of one of the altogether 6 mea¬ 
sured transmission spectra are shown in Fig. [T^ 



Frequency (GHz) 


FIG. 13: Transmission spectrum of the three-dimensional res¬ 
onator shown in Fig. for two different frequency ranges. 
Due to the high quality factors the resonances were well iso¬ 
lated up to 20 GHz yielding a complete sequence of altogether 
18764 resonance frequencies. 

For the investigation of the fluctuation properties of 
the resonance frequencies these had to be unfolded with 


Weyl’s formula for the smooth part NgmoothU) of the in¬ 
tegrated resonance density N{f) of electromagnetic cav¬ 
ities [59j |60]. It is a polynomial of third order in the 
frequency, where, in distinction to that for the corre¬ 
sponding quantum billiard, the quadratic term is absent. 
The fluctuating part, Nfiudf) = ^(/) “ ^smoothif)^ is 
shown in Fig.[^ Like in the case of the two-dimensional 
stadium billiard in Fig. it exhibits slow oscillations, 
that again are due to non-generic POs. These consist of 
marginally stable bouncing-ball orbits with lengths 2ri 
and 2 r 2 that bounce back and forth between the top and 
bottom walls of the quarter cylinders. Furthermore there 
are orbits in their common rectangular plane that are lin¬ 
early stable within the plane and unstable with respect 
to perpendicular disturbances. An analytical expression 
for the contribution of the former to the resonance den¬ 
sity was derived in [56], that for the entire non-generic 
one in [58] on the basis of a quantum adiabatic method. 
The latter is plotted as red (gray) line in Fig. 
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FIG. 14: (Golor online) Fluctuating part of the integrated re¬ 
sonance density (black line) obtained for the resonator shown 
in Fig. ^1 compared to the contribution of the non-generic 
orbits (red [gray] line). Reprinted from Phys. Rev. Lett. 89, 
064101 (2002). 


Figure \TE\ exhibits the experimental length spectrum, 
deduced from the Fourier transform of the fluctuating 
part of the density of the resonance frequencies as black 
full line. The red dashed one was obtained from the 
sum of the Fourier transforms of the analytical results 
for the resonance density of the non-generic POs and the 
unstable ones. The latter is provided by the semiclassical 
trace formula of Balian and Duplantier m, 


Pfluc 




Ltp 2 cos 

TT ^ |det(l - Mp)|i/2 


(j^Lp 2 ^p). 


(4) 

The factor 2cos(/)p accounts for the polarization, that 
is, the vectorial character of the Helmholtz equation for 
three-dimensional resonators. It vanishes for orbits with 
an odd number of reflections. Note, that for separa¬ 
ble systems, where the transversal electric and transver¬ 
sal magnetic modes decouple completely, |2cos0p| = 2. 
The quantities Lp and Mp denote the length of the peri- 
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odic orbit p, and its stability matrix, respectively. They 
are identical with those entering Gutzwiller’s trace for¬ 
mula [50]. The diamonds on the abscissa in Fig. 15 mark 
the peaks that are exclucively assigned to POs in the 
common plane of the quarter cylinders. The other peaks 
below / ^ 0.6 m are at the lengths of the bouncing-ball 
orbits. These indeed are comparatively large. 

In conclusion, the contributions from non-generic or¬ 
bits could not be neglected and had to be extracted for 
the study of the spectral fluctuation properties of the 
unfolded resonance frequencies. These were expected to 



FIG. 15: (Color online) Fourier transforms of the fluctuating 
parts of the r eson ance densities obtained for the resonator 
shown in Fig. 12 (black full line) and for the non-generic 
modes (58] |6T] plus the trace formula Eq. ffl of Balian and 


Duplantier (red [gray] dashed line). Both curves exhibit peaks 
at the lengths of the POs. The diamonds mark the lengths 
of the POs in the common rectangular plane of the quarter 
cylinders. 


coincide with the RMT result for the GOE, since the dy¬ 
namics of the corresponding classical billiard was shown 
to be chaotic [62[ [63]. In Fig. the experimental re¬ 
sult for the As statistics is plotted (diamonds) for three 
frequency ranges. Even though the non-generic contri¬ 
butions were extracted by proceeding as in the case of 
the two-dimensional stadium billiard (see Sec. Ell , we 
observe deviations from the GOE result (dotted lines). 
Instead, the As statistics is very well decribed by that 
of random matrices of the form in Eq. (§, applicable to 
sytems consisting of two coupled chaotic ones (full lines). 
Actually, in the frequency range of 2 — 5 GHz it is close 
to that of two uncoupled GOEs (dashed line). 

The reason for the deviations from the GOE results 
is not yet fully understood. Numerical studies (61] of 
the spectral properties of three-dimensional quantum sta¬ 
dium billiards demonstrated that they cannot be at¬ 
tributed to the partial breaking of the reflection symme¬ 
try present, when the radii ri and r 2 of the two quarter 
cylinders are equal. Furthermore, our interpretation of 
the good agreement of the spectral properties with those 
of the model in Eq. as a partial decoupling of the 
electric and the magnetic field modes could not yet be 
confirmed, neither experimentally nor numerically. How¬ 
ever, the fact that a considerable fraction of the POs have 


\2 cos (j)p\ = 2 provides a supportive indicator for its cor¬ 
rectness. In the following section we will present another 
system whose spectral properties deviate from the RMT 
predictions. 



EIG. 16: The A 3 statistics of the resonance frequencies of the 
resonator shown in Eig. [T^ evaluated in the three frequency 
intervals indicated in the insets (diamonds). The dotted curve 
corresponds to that for one GOE, the dashed one to that of 
two uncoupled GOEs, i.e., to the case A = 0 in Eq. (|^, and 
the full one to the curve deduced from Eq. § best fitting the 
experimental one yielding A = 0.005, 0.124, 0.498 from top to 
bottom. 


VI. SPECTRAL FLUCTUATIONS OF A 
SYSTEM WITH THREEFOLD SYMMETRY 

As has been already stated above, according to 
the Bohigas-Giannoni-Schmit conjecture [5| the spectral 
properties of time-reversal (T) invariant chaotic systems 
coincide with those of random matrices drawn from the 
GOE. This, however, is not true for systems with an uni¬ 
directional classical dynamics (see Sec. Ivnl) or for bil¬ 
liards that have a shape with a threefold symmetry and 
no further one like, e.g., a mirror symmetry |64|. In this 
context, we investigated the spectral properties of the 
microwave billiard presented in Fig. Its shape is in¬ 
variant under rotations of 120°. The eigenfunctions T of 
the Hamiltonian of such a quantum billiard can be classi¬ 
fied according to their transformation properties with re¬ 
spect to a rotation R by 120°, RT = exp (27ri//3) T, I = 
— 1,0,1. Due to Tinvariance the eigenvalues associated 
with the classes I = ±1 are degenerate. Accordingly, 
the eigenvalue spectrum is composed of singlets corre¬ 
sponding to the case I = 0 and doublets of identical 
ones [64l [65] . The spectral properties of the singlets and 
the doublets have been predicted to coincide with those 
of random matrices from the GOE and the GUE, respec¬ 
tively [64]. 
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FIG. 17: (Color online) Photograph of the microwave billiard 
with threefold symmetry and no further symmetries. It con¬ 
sisted of a bottom plate, a middle plate, which had a hole 
with the shape of the billiard, and a lid. All parts were con¬ 
structed from copper and covered with lead before they were 
squeezed together tightly with screws along the boundary of 
the cavity. The top plate is removed. 


We determined the eigenvalues of the billiard with the 
shape shown in Fig. [^by performing measurements at 
superconducting conditions. Still, since the degeneracy 
of the doublet states was only partly split due to mechan¬ 
ical imperfection [ 66 ] a more sophisticated method had 
to be developed for their identification m- It utilized 
the dependence of the behavior of the eigenmodes un¬ 
der the variation of the relative phase between two input 
antennas on their symmetry properties. In Fig. [Ts] the 
nearest-neighbor spacing distributions of the singlets and 
doublets (histograms) are compared with the RMT result 
for the GOE (black line) and for the GUE (red [gray] 
line), respectively. Similarly, the A 3 statistics of the sin¬ 
glets and doublets (diamonds) is presented in Eig. [l^ 
together with the GOE (black line) and the GUE (red 
[gray] line), respectively. These figures demonstrate that 
the experimental results were in accordance with the the¬ 
ory. 


VII. DYNAMICAL TUNNELING IN 
CONSTANT-WIDTH BILLIARDS 

Recently, we investigated the phenomenon of dynami¬ 
cal tunneling in two microwave billiards with the shapes 
of constant-width billiards [ 68 l [69| . A photograph of one 
is shown in Eig. Such billiards possess the particular 
property of unidirectionality, that is, when a particle is 
launched into the billiard with a certain rotational direc¬ 
tion then it will continue to rotate in that direction for¬ 
ever. Accordingly, when constructing a Poincare surface 
of section (PSOS) by launching many particles into the 
billiard in clockwise direction, then only its upper half 
will be filled, as illustrated in Eig. Here, p = sin 6 > 
where 0 is the angle between the particle trajectory and 
the normal to the boundary at the point of impact, and 
q is the location of the latter in terms of the arclength 



FIG. 18: (Golor online) Nearest-neighbor spacing distribu¬ 
tions of the resonance frequencies of the singlets (upper panel) 
and the doublets (lower panel) obtained for the microwave bil¬ 
liard shown in Fig. [Ff| (histogram). The RMT results for the 
GOE and the GUE are shown as black and red (gray) lines, 
respectively. 



FIG. 19: (Golor online) The A 3 statistics of the resonance 
frequencies of the singlets and the doublets obtained for the 
microwave billiard shown in Fig. [T7| (diamonds). The RMT 
results for the GOE and the GUE are plotted as black and 
as red (gray) lines, respectively. Reprinted from Phys. Rev. 
Lett. 90, 014102 (2003). 

along the boundary. The PSOS consists of a large chaotic 
sea and two regions of regular KAM tori. Eirstly, there 
is a barrier region around the diameter orbit at p = 0 , 
shown in blue (light gray). Secondly, there is a region 
of whispering gallery modes around p = 1. Examples 
of both types of orbits are plotted in the inset in the 
corresponding colors. Eurthermore, we found very tiny 
regular islands inside the chaotic sea for one of the two 
investigated billiards m- 
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FIG. 20: (Color online) Photograph of the microwave billiard 
with the shape of the constant-width billiard Bl. It was con¬ 
structed from copper plates. The cavity was milled out of 
the bottom plate. The top plate is removed. Reprinted from 
Phys. Rev. E 90, 022903 (2014). 


widths of the barrier and the whispering gallery orbit re¬ 
gions were considerably larger in the PSOS of Bl than in 
that of B2. More details concerning the measurements 
are provided in Ref. [19]. 

To demonstrate the high precision of the measurements 
and also its indispensability for the identification of the 
doublets we show in Fig. two parts of the transmission 
spectrum of billiard Bl. The resonances visible in the 



The fact, that only one half of the PSOS is filled is a 
clear indication that classically the transition from clock¬ 
wise to anticlockwise motion, i.e., the transition from 
positive to negative p values is strightly forbidden. In 
contrast, in the corresponding quantum billiard this is 
possible via tunneling through the barrier region around 
p = 0. This manifests itself in a partial lifting of the de¬ 
generacies of pairs of eigenenergies associated with the 
clockwise and anticlockwise quasimodes [T9|. Accord¬ 
ingly, the size of the tunneling effect is reflected in that 
of the splittings of the eigenenergies. 


ra 

II 


FIG. 21: (Golor online) Poincare surface of section (PSOS) for 
billiard Bl where all particles were launched into the billiard 
in the clockwise direction. Due to unidirectionality only the 
upper half of the PSOS is filled. The chaotic sea is bordered 
by a region of whispering gallery modes around p = 1 plotted 
in red (gray) and one of orbits around the diameter orbit 
plotted in blue (light gray). The inset shows examples for the 
associated orbits. Reprinted from Phys. Rev. E 90, 022903 
(2014). 



The objective of the experiments with the two mi¬ 
crowave billiards with the shapes of constant-width bil¬ 
liards, referred to as billiards Bl amd B2 in the se¬ 
quel, was to investigate this phenomenon with unprece¬ 
dented accuracy using superconducting microwave bil¬ 
liards. Both billiards had the same diameter 2r = 48 cm 
and their shapes were barely distinguishable, whereas the 


-20 
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Frequency (GHz) 


EIG. 22: Resonance spectrum of billiard Bl. The upper panel 
shows it in a frequency range from 15.4 to 15.9 GHz, the 
lower one a zoom into a triplet of nearly degenerate pairs of 
resonances. Reprinted from Phys. Rev. E 90, 022903 (2014). 


upper panel seem to be sharp and well isolated. However, 
a zoom into the region of the triplett observed at the left 
side of the spectrum, reveals that they all correspond 
to doublets of nearly degenerate resonance frequencies. 
In total, 390 doublets and 15 singlets were detected in 
billiards Bl and B2, respectively. The latter correspond 
to diameter orbits, that are analogues of those of the 
circle billiard of the same diameter 2r, the eigenvalues 
are the zeroes of the Bessel functions, Jo{k^r) = 0. 
We actually identified all singlets in the list of resonance 
frequencies of Bl and B2 by comparison with these zero- 
momentum modes. 

In order to check the completeness of the sequence 
of doublets and to study the spectral properties of the 
constant-width billiards, we separated the resonance fre¬ 
quencies into singlets, and the smaller and larger ones 
of the doublets, i.e., the left, /^, and the right, /^, dou¬ 
blet partners in the resonance spectra (see lower panel of 
Fig. 22), respectively. This was possible, because the cor¬ 
relations between the doublet partners and those between 
the singlets and the doublets were much weaker than 
those amongst the members of each of the sequences. 
Consequently, we could simply disregard the singlets and 
investigate the fluctuation properties of the left and the 
right doublet partners separately. For this, we first un¬ 
folded them with Weyl’s formula [21] (see Sec. IIIip . The 
comparison of the smooth part of the integrated density 
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of the resonance frequencies with the latter and also the 
shape of its fluctuating part confirmed that we identified 
the complete sequence of doublets. Due to the unidi¬ 
rectionality of the classical dynamics, the spectral prop¬ 
erties were predicted to coincide with the RMT result 
for the GUE, even though the system is Tinvariant [69] . 
We didn’t find any agreement with the GUE and also not 
with the GOE results, but instead found very good agree¬ 
ment with those of a random-matrix model for systems 
with violated Tinvariance and mixed dynamics. Within 
this model the form of the random matrices is given by 
Eq. (§, with Hi replaced by a random matrix from the 
GUE, and H 2 by a diagonal matrix with random Pois- 
sonian numbers as entries. The elements of the matrix 
U, which describes the mixing of the chaotic and the 
regular states modeled by Hi and i72, respectively, were 
chosen as Gaussian-distributed random numbers. Thus, 
even though the classical dynamics of the investigated 
billiards was not fully chaotic, the predictions of Ref. [69] 
could be confirmed. 

— I — I — I — I — I —I— I —I— I — 



Frequency (GHz) 

FIG. 23: (Color online) The splittings of the resonance fre¬ 
quencies of the doublets, S = /^ —//^, for B1 (upper panel) and 
B2 (lower panel) versus their average resonance frequencies. 
They are by a factor of 10^ — 10^ smaller than their average 
resonance frequency. This underlines again the necessity of 
high-precision experiments to identify a complete sequence of 
doublets. 

The splittings S = — of the resonance frequencies 

of the doublet partners are plotted versus their average 
resonance frequency for B1 in the upper panel and for 
B2 in the lower one of Eig. respectively. The typical 
splitting size is smaller for B1 than for B2. This is in 
agreement with the expectation, that the tunneling rate 
through the barrier region in the PSOS, which is broader 
for B1 than for B2, decreases with increasing width of 
that region. Eurthermore, the average splitting size di¬ 
minishes algebraically with increasing frequency, not ex¬ 
ponentially as might be expected semiclassically. This 
can be attributed to the fact, that the PSOS consists 
of a chaotic region and several regular regions, where the 


splitting sizes diminish exponentially with different rates. 
We considered all splittings, i.e., the sum over the expo¬ 
nentials which yields an algebraical decay 170]. 




FIG. 24: (Color online) Comparison of the experimentally 
determined splitting distributions (black histograms) of B1 
[upper panel] and B2 [lower panel] with that obtained from 
the random-matrix model of Eq. ra shown as red (gray) his¬ 
tograms and with the analytical result of Eq. 0 plotted as 
green (light gray) line. Here, the tunneling parameters de¬ 
duced from the fit of the RMT model to the experimental 
distributions are r = 0.094 for B1 and r = 0.185 for B2. 
Reprinted from Phys. Rev. E 90, 022903 (2014). 

Eor a better understanding of the tunneling we inves¬ 
tigated also the distribution of the splittings. The black 
histograms in the upper and lower panel in Eig. [^corre¬ 
spond to the experimental results for billiards B1 and B2, 
respectively. In order to extract system specific quanti¬ 
ties, like the area and the perimeter of the billiards, we 
used the splittings S of the unfolded left and right doublet 
partners. We demonstrated in Ref. m that the shapes 
of the splitting distributions are not affected by states 
that are associated with the regular parts of the PSOS, 
that is, the chaotic component dominates. They, indeed, 
are well described by a random-matrix model of the form 
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which is applicable to systems with broken Tinvariance 
exhibiting tunneling between two regions with chaotic 
dynamics. Here, H and Ff* are 2x2 random matrices 
from the GUE, that are associated with the parts of the 
Hilbert space corresponding to clockwise and anticlock¬ 
wise motion, respectively, and Dq is a random Poissonian 
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number which represents a singlet state. The tunneling 
from the upper to the lower part of the PSOS via the 
barrier region, that is, the coupling of H and i7* via 
the singlet state, is caused by the off-diagonal matrix V, 
which has Gaussian-distributed random numbers as en¬ 
tries. The strength of the coupling is determined by the 
tunneling parameter r. It was obtained by fitting the 
splitting distribution of the RMT model of Eq. (§ to the 
experimental ones. On the basis of the latter we further¬ 
more derived an analytical expression for the splitting 
distribution. 


p(5) = 


2 


-4 

e 


r4 

(5 + r2)2 


1 + 8 


((5 + r2)2j 


, ( 6 ) 


which we believe is generally applicable to systems with 
violated Tinvariance, where dynamical tunneling takes 
place between two regions with a chaotic dynamics. The 
green (light gray) lines in Fig. 24 correspond to the an- 
alyticial results computed for the parameter values of r 
given in the captions. Its value is smaller for billiard B1 
than for B2, as expected, because the PSOS of the for¬ 
mer comprises a broader barrier region. The agreement 
between the experimental, the RMT and the analytical 
results for the splitting distributions is very good. Re¬ 
markably, the RMT model and the analytical formula for 
the splitting distribution depend only on one parameter 
which provides information on the tunneling rate. 

We would like to close this section by mentioning an¬ 
other experiment on dynamical tunneling, which we per¬ 
formed 15 years ago with superconducting microwave 
billiards with the shapes of annular billiards [HI E]. 
There, the dynamical tunneling between quasimodes cor¬ 
responding to clockwise and anticlockwise moving whis¬ 
pering gallery orbits, respectively, was investigated. As 
recognizable in the upper panels of Figs. and such 
a billiard corresponds to a circular one of radius R that 
contains a circular disk of radius r < R, where the cir¬ 
cle centers are displaced by a value S with respect to 
each other. The classical dynamics is characterized by 
the value (r + S)/R which was chosen equal to 0.75. 
Then the regions of neutrally stable clockwise and an¬ 
ticlockwise moving whispering gallery orbits are located 
in the PSOS at angular momentum values S > 0.75 and 
S < —0.75 in units of R, respectively, where S = sinO 
corresponds to p in Fig. 21 For |5'| = | sin^| < 0.75 the 


PSOS consists of a large chaotic sea with a regular island 
at its center. 

The objective of the experiments was the investigation 
of chaos-assisted tunneling [721173] between quasimodes 
localized in the regions |5'| > 0.75 via the chaotic sea, 
that is, the enhancement of the tunneling through a cou¬ 
pling of these modes to chaotic ones. The coupling is 
largest for modes located close to the borders S ^ ±0.75 
that separate the three regions. As in the case of the 
constant-width billiard the tunneling manifests itself in 
the splitting of the degenerate quasimodes corresponding 
to clockwise and anticlockwise moving whispering gallery 
orbits into doublets. Thus, if chaos-assisted tunneling 


takes place, its size should depend on their location in 
the PSOS, i.e., on their angular momentum values. Ac¬ 
cording to the Einstein-Brillouin-Keller quantization the 
quantum analogue of the angular momentum S = sin^ 
is given as S' = n/kn,m for a pair of whispering gallery 
modes. Here, kn,m denotes their mean wave number and 
their angular and radial quantum numbers are n and m, 
respectively. They actually coincide to a good approxi¬ 
mation with those of the integrable annular billiard with 
(5 = 0. To determine them the electric field distributions, 
i.e. the wave functions of the annular billiard, were mea¬ 
sured in a normal conducting microwave billiard made 
from copper using the perturbation method [niEl]. The 
mean wave numbers kn,m of the whispering gallery modes 
and the associated splittings were obtained from reso¬ 
nance spectra that were taken with superconducting bil¬ 
liards constructed from niobium. Figureshows a small 
part of a transmission spectrum which contains two sin¬ 
glets and one doublet, which is assigned to a pair of whis¬ 
pering gallery modes with the angular quantum number 
n = 18 given by half the number of maxima of the field 
intensity shown in the middle panel above the spectra. 
Furthermore, it has radial quantum numer m = 1, since 
the latter exhibits only one maximum in the radial direc¬ 
tion. The singlets correspond to modes that are localized 
in the region with |5'| < 0.75. 



FIG. 25: (Color online) Transmission spectrum of the annu¬ 
lar billiard shown schematically in the three top panels in 
the frequency range from 8.88 to 8.96 GHz obtained at room 
temperature (upper panel) and at 4.2 K (lower panel). The 
combined measurement of field distributions with the normal 
conducting microwave billiard and of resonance spectra with 
the superconducting one allows the identification of the quasi¬ 
doublets and the assigment of quantum numbers (m, n) and 
of a parity with respect to the symmetry axis like for the 
three modes associated with the resonances, shown in the top 
panels. Reprinted from Phys. Rev. Lett. 84, 867 (2000). 


For the study of the dependence of the tunneling rates, 
that is, of the sizes of the splittings on the location of 
the associated doublet modes in the PSOS, i.e., on their 
angular momentum values, we considered only modes 
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with radial quantum number m = 1. Figure shows 
the resulting splittings versus the angular momentum 
S = n/k^^ with = 27r/^/c denoting the wave num¬ 
ber corresponding to the resonance frequency of the 
mode. The curve exhibits a clear peak in the splittings, 
i.e., an enhancement of the tunneling in the beach region 
around the border |5'| = 0.75 between the regular coast 
and the chaotic sea. To the best of our knowledge, this 
corroborated for the first time the occurrence of chaos- 
assisted tunneling. 



FIG. 26: (Color online) Splittings of the doublets of reso¬ 
nances (squares) with respect to their average resonance fre¬ 
quency versus its position in the PSOS. The error in the 
determination of the frequencies of the doublet partners is 
generally smaller than the size of the squares except for the 
smallest splitting observed. A smooth transition from large 
splittings in the chaotic sea {S < 0.75) towards very small 
ones in the regular coastal region {S > 0.75) is visible. The 
appearance of a pronounced maximum in the direct vicinity 
of the beach {S ~ 0.75) is a clear indicator for the occurrence 
of chaos-assisted tunneling. Reprinted from Phys. Rev. Lett. 
84, 867 (2000). 


VIII. SPECTRAL PROPERTIES OF 
SUPERCONDUCTING MICROWAVE DIRAC 
BILLIARDS MODELLING GRAPHENE 

Only recently we have started a series of experiments 
with microwave photonic crystals [75] that model prop¬ 
erties of graphene ngizzi, a monoatomic layer of carbon 
atoms arranged on a honeycomb lattice. Due to its par¬ 
ticular electronic properties, graphene attracted a lot of 
attention over the last years. Its conductance and its va¬ 
lence band form conically shaped valleys that touch each 
other at the six corners of the hexagonal Brillouin zone. 
In their vicinity the electron energy depends linearly on 
the quasimomentum, implying an energy independent ve¬ 
locity. Thus, the electrons and holes behave like massless 


relativistic particles governed by the Dirac equation for 
spin-1/2 particles. Accordingly, the touching points are 
referred to as Dirac points [iniiii]. 

This extraordinary feature of the band structure of 
graphene results in a number of electronic properties 
which have an analogue in relativistic quantum mecha¬ 
nics. It actually stems from the symmetry properties 
of its honeycomb structure which is formed by two in¬ 
terpenetrating triangular lattices with threefold symme¬ 
try. This led to the manufacturing of artificial graphene 
using two-dimensional electron gases, molecular assem¬ 
blies, ultracold atoms [784182] or photonic crystals [Sl- 
ED]- We exploited this fact to simulate the spectral prop¬ 
erties of graphene with microwave resonators that were 
constructed by squeezing a photonic crystal, which is 
composed of several hundreds of metallic cylinders ar¬ 
ranged on a triangular lattice, between two metal plates. 
The shape of its band structure, i.e., of its frequencies of 
wave propagation as function of the two components of 
the quasimomentum vector, is equivalent to that of the 
conductance and the valence band in graphene, e.g., for 
its 1st and its 2nd band. Accordingly, around the Dirac 
points the propagation of the microwaves in the photonic 
crystal can be effectively described by the Dirac equation 
for massless spin-1/2 particles. Such a region is identified 
in the resonance spectrum as one of an exceptionally low 
resonance density, which increases linearly with the dis¬ 
tance from the frequency of the Dirac point, the so-called 
Dirac frequency {83] [84] . 

In order to investigate this behavior in detail, we per¬ 
formed high-precision experiments with superconducting 
Dirac billiards with the shapes of a rectangle and the 
African continent [9T4I93] . respectively. Photographs are 
shown in the left parts of Figs. and [^ We use them 
to simulate properties of graphene billiards or graphene 
flakes. Both billiards consist of a bassin with the shape 
of the billiard and a lid, that were made from brass plates 
and then lead plated to attain superconductivity at liquid 
helium temperature. For the construction of the photonic 
crystals located inside the bassins ~ 900 cylinders were 
milled out of a brass plate. The height of the resonators 
was d = 3 mm so for / < fmax = 50 GHz the associated 
Helmholtz equation is mathematically identical to the 
Schrodinger equation of the quantum billiard containing 
circular scatterers at the positions of the cylinders. 

The great advantage of Dirac billiards as compared 
to graphene flakes is that we can control the boundary 
conditions. Here, those for the wave functions at the 
voids that are located at the edges of the photonic crystal 
are of relevance. They can be manipulated by varying 
the relative position between these voids and the walls of 
the microwave billiard, where the electric field strength 
vanishes, i.e., obeys the Dirichlet boundary condition. 

The cylinders were arranged on a triangular grid, as in¬ 
dicated in the schematic views in Figs. and so the 
voids at the centers of the cells formed by, respectively, 
three of them yield a hexagonal configuration similar to 
the honeycomb structure in graphene. Each of the cells 
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FIG. 27: (Color online) Left Part: Photograph of the Dirac 
billiard with the shape of a rectangle, which contains 888 
metal cylinders arranged on a triangular grid. The top plate 
is shifted with respect to the bottom one. Reprinted from 
Phys. Rev. B 88, 104101 (2013). Right Part: Schematic view 
of the triangular lattice of the rectangular Dirac billiard. The 
red (gray) and blue (dark gray) dots mark the voids between 
the cylinders. They form a honeycomb lattice. 



FIG. 28: (Color online) Same as in Fig. 27 but for an Africa 
billiard. The top plate is removed. 


may host a quasibound state and thus can be considered 
as one of the atoms forming the two independent triangu¬ 
lar sublattices in graphene. Accordingly, both graphene 
and microwave photonic crystals can be treated in a sim¬ 
ilar way within the tight-binding model (TBM). 

We determined altogether 1651 and 1823 resonance fre¬ 
quencies in the first two bands of the rectangular and 
the Africa Dirac billiard, respectively. Figure shows 
their resonance densities (black curves). They clearly 
differ from those expected for the corresponding empty 
microwave billiards, that, according to WeyFs law, would 
increase linearly with the resonance frequency. Both ex¬ 
hibit a minimum at the frequency of the Dirac point, 
which separates the two bands. It is bordered by two 
sharp peaks corresponding to van Hove singularities, that 
generally occur in two-dimensional crystals with a peri¬ 
odic structure [95]. There, the resonance density diverges 
logarithmically for infintely extended crystals, whereas in 
the case of the Dirac billiards the peak heights increase 
logarithmically with the size of the crystal, i.e. with the 
number of voids formed by the cylinders [96] . 


The red (gray) dashed lines in both panels of Fig. 29 


correspond to fits of the TBM m for finite-size pho¬ 
tonic crystals to the experimental ones. It incorporates 
nearest-neighbor, as well as second- and third-nearest- 
neighbor couplings and in addition the corresponding 
overlaps of the wave functions centered at the different 
voids. Furthermore, the resonance density of the Africa 
Dirac billiard exhibits a slight bump to the right of the 


Dirac frequency, which is attributed to the zigzag edge 
parts of the void structure formed by the associated pho¬ 
tonic crystal [94l 198] [99| . The voids of the rectangular 
billiard also have a zigzag structure along its longer sides, 
see Fig. ^7\ however, the photonic crystal can be repro¬ 
duced periodically via reflections at all its sides so no 
edge states are generated [94|. 



FIG. 29: (Golor online) Densities of the resonance frequencies 
of the Dirac billiards with the shape of a rectangle [upper 
panel] and of the African continent [lower panel]. A slight 
bump is visible in the latter to the right of the Dirac frequency. 
It is attributed to edge states occurring along the zigzag edge 
parts of the void structure, see also the right hand side of 
Fig.m 


We demonstrated recently in Ref. [94| that the van 
Hove singularities divide the two bands framing the Dirac 
point into two frequency ranges below and above them, 
where the system is described by the nonrelativistic 
Schrodinger equation of the empty quantum billiard, and 
one between them, where the underlying wave equation 
coincides with the Dirac equation of the graphene bil¬ 
liard or flake of corresponding shape. Actually, the TBM 
does not only provide the eigenvalues, but also the wave 
functions of a bounded photonic crystal. In the regions 
close to the lower and the upper band edges, they are 
identieal with those of the corresponding quantum bil¬ 
liard. Moreover, we showed that this also holds for the 
fluctuating part of the integrated resonance density and 
for the spectral fluctuation properties. Furthermore, the 
computation of the momentum distributions using the 
wavefunctions and also the direct comparison of the band 
structure of the photonic crystals with the resonance fre¬ 
quencies plotted versus the eigenvalues of the quantum 
billiard yielded, that in the vicinity of the lower and the 
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upper band edge the latter maybe identified with the 
quasimomenta. Accordingly, we call the regions below 
and above the van Hove singularities the nonrelativistic 
Schrodinger region. 

In the region between the van Hove singularities, i.e., 
around the Dirac point, the electromagnetic waves inside 
the billiards are governed by the Dirac equation. Ac¬ 
cordingly, there we compared the resonance frequencies 
and their fluctuation properties with those of the eigen¬ 
values of the corresponding graphene billiard [98]. The 
latter were computed by solving the Dirac equation with 
appropriate boundary conditions depending on the edge 
structure. We arrived at the conclusion, that they and 
the resonance frequencies interdepend linearly and fur¬ 
thermore, that they coincide with the quasimomenta in 
the vicinity of the Dirac point. Thus, at the van Hove sin¬ 
gularities a transition takes place between a nonrelativis¬ 
tic Schrodinger region and the relativistic Dirac region, 
where the Dirac billiard exhibits the same features as the 
corresponding empty quantum billiard and the graphene 
billiard, respectively. In Ref. [96] we provided evidence 
that it can be interpreted as a Lifshitz and an excited 
state quantum phase transition. Presently, we perform 
more investigations concerning the spectral properties of 
the Dirac billiards at the van Hove singularities in order 
to corroborate this assumption. 

The spectral properties are similar in the different re¬ 
gions and agree with those of the corresponding empty 
quantum billiard. The classical dynamics of the empty 
rectangular and Africa billiards is regular, respectively, 
chaotic [911493] . The nearest-neighbor spacing distribu¬ 
tion and the A 3 statistics of the rectangular Dirac bil¬ 
liard shown in the upper panels of Fig. (histogram) 
and Fig. (diamonds), respectively, are well described 
by those of Poissonian random numbers (red [gray] lines) 
whereas those of the Africa one plotted in the lower pan¬ 
els coincide with the GOE results (black line). 


IX. CONCLUSIONS 

The objective of this article was to present some of the 
highlights of our experimental investigations concerning 
the quantum manifestations of classical chaos. The focus 
has been set on spectral fluctuation properties, on sym¬ 
metry breaking showing up in particular properties of the 
eigenvalues and the wave functions, on the phenomenon 
of dynamical tunneling, and, moreover, the universal 
properties of graphene. Due to the limited space, we ref- 
ered to only a few examples of the numerous experiments 
performed with superconducting microwave cavities in 
the last two decades. Accordingly, we selected some that 
provided groundbreaking new insight into physical phe¬ 
nomena. Future projects concern the realization of scat¬ 
tering systems by using microwave billiards with one or 
several openings in the side walls [10011102] . The objec¬ 
tive is the measurement of the temporal decay of mi¬ 
crowave power through these openings and a comparison 



FIG. 30: (Color online) Nearest-neighbor spacing distribu¬ 
tions of the resonance frequencies obtained for the Dirac bil¬ 
liards (histograms) with the shapes of a rectangle (upper 
panel) and of the African continent (lower panel). The former 
agrees well with the Poisson distribution (red [gray] line), the 
latter with the RMT result for the GOE (black line). The up¬ 
per figure has been reprinted from Phys. Rev. B 91, 035411 
(2015). 



FIG. 31: (Color online) The A 3 statistics of the experimen¬ 
tally determined resonance frequencies of the Dirac billiards 
(diamonds) with the shapes of a rectangle (upper panel) and 
of the African continent (lower panel). The former agrees 
well with that of Poissononian random numbers, the latter 
with the RMT result for the GOE. The upper figure has been 
reprinted Phys. Rev. B 91, 035411 (2015). 


with RMT and other predictions. Such experiments yield 
information on the universal decay behavior of open sys¬ 
tems like excited nuclei. In order to ensure that the sys¬ 
tem decays exclusively via the openings, the dissipation 
into the walls of the resonator needs to be eliminated. 
This is only possible by performing the measurements 















17 


of the S'-matrix elements at superconducting conditions. 
There, however, a calibrated measurement, an essential 
requirement for the analysis of properties of the S matrix, 
is not yet feasible. Such a method was developed only 
recently Currently we implement and adjust it to 

our superconducting measurements of the S matrix with 
open microwave resonators with the shapes of regular 
and chaotic billiards. The aim is to obtain information 
on the differences in their decay behavior. Furthermore, 
we construct superconducting quantum graphs [ina, a 
network-shaped structure of vertices bonded by rectan¬ 
gular waveguides. The objective of these experiments 
is the experimental test of RMT predictions concerning 
the correlation functions of the S matrix of open quan¬ 
tum graphs |1Q5[ 1106] in high-precision experiments. Fi¬ 
nally, a further current project concerns the experimen¬ 
tal verification of the so-called Atiyah-Singer index the¬ 
orem fTTmmn] . which relates the number of zero modes 
of the Dirac operator to the topology of the associated 
graphene sheet. The resonator consists of a network of 
rectangular waveguides that connect circular cavities. It 


was milled out of a brass sphere and has the form of a 
fullerene Ceo molecule m- Summarizing, we have with 
our macroscopic, superconducting microwave resonators 
systems available that might be tailored to model the uni¬ 
versal properties of a large range of microscopic quantum 
systems in high precision experiments. 
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